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ABSTRACT 



We present initial results of an ESO-VLT large programme (AMAZE) aimed at determining the evolution of the mass-metallicity 
relation at z>3 by means of deep near-IR spectroscopy. Gas metallicities are measured, for an initial sample of nine star forming 
galaxies at z~3.5, by means of optical nebular lines redshifted into the near-IR. Stellar masses are accurately determined by using 
Spitzer-IRAC data, which sample the rest-frame near-IR stellar light in these distant galaxies. When compared with previous surveys, 
the mass-metallicity relation inferred at z~3.5 shows an evolution much stronger than observed at lower redshifts. The evolution is 
prominent even in massive galaxies, indicating that z~3 is an epoch of major action in terms of star formation and metal enrichment 
also for massive systems. There are also indications that the metallicity evolution of low mass galaxies is stronger relative to high mass 
systems, an effect which can be considered the chemical version of the galaxy downsizing. The mass-metallicity relation observed at 
z~3.5 is difficult to reconcile with the predictions of some hierarchical evolutionary models. Such discrepancies suggest that at z>3 
galaxies are assembled mostly with relatively un-evolved sub-units, i.e. small galaxies with low star formation efficiency. The bulk of 
the star formation and metallicity evolution probably occurs once small galaxies are already assembled into bigger systems. 

Key words. ISM: abundances - galaxies: abundances - galaxies: evolution - galaxies: high-redshift - galaxies: starburst 



> 1. Introduction 



The connection between galaxy mass and metallicity has been 
known for a long time, starting with the seminal work of 
iLequeux et alj dl979l) . Given the difficulty in obtaining reliable 
galaxy masses several authors have resorted in using the (opti- 
cal) luminosity. In particular, various works have reported a clear 
correlation between blue luminosity and metallicity, in the sense 
that more lumi nous galaxies are charac t erized by higher metal- 
licities (e.g. iGarnett & Shieldsl [19871: ISkillman et all Il989t 
iBrodie & Huchralll99ll:IZaritskv et allll994l). A m a jorste p for- 
ward has been recently achieved bv lTremonti et all 0004). who 
used optical photometric and spectroscopic data of a sample 
of ~53,000 galaxies from the SDSS to determine the mass- 
metallicity relation of local galaxies (z~0.1). Their work clearly 
showed that the primary physical parameter driving the corre- 
lation with the gas metallicity i s the (stellar) mass o f galaxies 
and not their luminosity. While Tre monti et alj d2004l) focused 
on the gas metallicity, a similar relation by using the SDSS sur- 
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vey was found by Gallazzi et al. (2006) for what concerns the 
stellar metallicity. 

Various physical processes may be responsible for the mass- 
metallicity relation. One possibility is that outflows, generated 
by starburst winds, eject metal-enriched gas into the IGM pref- 
erentially out of low-mass galaxies (due to the shallow gravita- 
tional potential well), ma king their enrichment less effective than 
in massive systems (e.g. iTremonti et all l2004t |Pe Lucia et all 
12004 |Finlator & Pavel I2008I) . An alternative scenario is that 
low mass systems are still at an early evolutionary stage and 
have still to convert most of their gas into stars, hence they 
are poorly metal-enriched relative massive galaxies (which are 
instead already evolved). This is the so-called "galaxy down- 
sizin g" scenario, supported by various observational evidences 
(e.g.lJuneau et alll2 005; Feul ner et alll2005HFranceschini et all 
2006; lAsari et all 120071: iPerez-Gonzalez et all 120071) . where 
massive galaxies formed most of their stars rapidly and at high 
redshift, while low mass systems are characterize d by a slower 
evolut ion, which extends to low redshift. Finally Koppen et al. 
(120071) ascribes the mass-metallicity relation to variations of the 
IMF high-mass cutoff in different star forming environments. 

The relative role of these processes in shaping the mass- 
metallicity relation is debated. It is likely that each of them 
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contributes at least to some extent, since observational ev- 
idences have been found for all of them. Each of these 
factors (outflows/feedback, downsizing, IMF) has profound 
implications on the evolution of galaxies. Therefore, it is 
clear that the mass-metallicity relation contains a wealth 
of information useful to constrain models of galaxy forma- 
tion and evolution. Indeed, any model of galaxy evolution 
is now required t o match the mass-metallicity relation ob- 
served locally (e.g. Kobayas hi et all 120071 IBrooks et all [2 007; 
de Rossi et all 2007 : iPave & Oppenheimej. l2007t iPalcantonl 



2007): iDe Lucia et all 120041 iTissera et all 120051; iBouche et al 



2006. 120071; iKoppen et all 120071; ICid FernandesetaU [2007; 
Finlator & Pavel [2008; Tas sis et all 120081) . However, differ- 



ent models predict different evolutionary patterns of the mass- 
metallicity relation as a function of redshift, and observational 
data are required to test and discriminate among them. 

Observational constraints of the mass-metallicity relation 
have been obtained up to z~2.2 thanks to various deep sur - 
veys (ISavaglio et all 120051: iLiang et all 120061: lErb et all [2006). 
Additional observational studies have investigated the evolu- 
tion of the luminosity-metallicity relation or , more generally, the 
metall icity of high-z star formin g galaxies ( Kobulnickv & Koo, 
2000; Kobulni ckv etail [20031 : TKobulnic kv & Kewlevl 12004; 
Maier et all 12004 I2005L (2006t iForster Schreiber et all [2006). 
Note that all of these studies refer to the gas metallicity, 
while limited work has been done on the stellar metalhcity 
of high redshi f t sources dRix et all 12004 Ide Mello et all 12004 
iHalliday et all 120081: iMehlert et al.[ l2006) due to difficulties in 
obtaining high S/N spectra on the stellar continuum. The gen- 
eral observational result is that the mass-metallicity relation 
(as well as the luminosity-metallicity relation) evolves, in the 
sense that on average higher redshift galaxies are character- 
ized by lower metallicities (at a given mass). Whether the re- 
lation evolves also in terms of its shape is still matter of debate. 
Theoretical models can generally cope with the observed evolu- 
tion, within both the t heore t ical and observational uncertainties 
(e.g. [ Kobayas hi et all 12007]: I de Rossi et all 120071 : IBrooks et all 
120071: iFinlator & Pavel 12008b . 

At z>3, except for a few individual studies, little work has 
been c urrently done for w hat concerns the mass-metallicity re- 
lation. iPettini et all d200ll) have measured the metallicity for 
a small sample of Lyman Break Galaxies (LBG) at z~3, 
but without investigating the mass-metallicity relation. The 
metalhcity evolution h a s been investigated in PLA systems 
dProchaska et all 120031: iKulkarni etail 120051: lAkerman etail 
120051) . However, a study of the mass-metallicity relation for 
high-z absorption systems has not been properly performed, due 
to difficulties in measuring the associated stellar masses. The 
closest result is the finding of a relation between metallicity 
and velocity dispersion (probably related t o the mass) inferred 
by the width of the a bsorption systems dLedoux et all 2006; 
iProchaska etail 120071) . 

We have undertaken a large observing programme at ESO- 
VLT aimed at carefully determining the mass-metallicity rela- 
tion at z>3 for a sizeable sample of galaxies. The final goal is 
to obtain a detailed description of the evolution of the mass- 
metallicity evolution through the cosmic epochs, and therefore 
constrain galaxy evolutionary scenarios. In this paper we present 
preliminary results obtained by such a program, and discuss 
the relevant implications for our understanding of the galaxy 
evolution at high redshift. Throughout the paper we adopt the 
following cosmological p arameters: Hp = 71 km s _1 Mpc -1 , 
Q A = 0.73 and Q m = 0.27 dSpergel et alll2003l) . 



2. The AMAZE program 

AMAZE (Assessing the Mass-Abundance redshift[-Z] 
Evolution) is an ESO large program aimed at determin- 
ing the mass-metallicity relation in the redshift range 
3<z<5. Observations are being performed with SINFONI 
(Eisenhau er et alll2003l) . the near-IR integral field spectrometer 
at VLT, for a total of 180 hours, distributed in three semesters. 
Observations are expected to be completed in mid-2008. The 
target sample consists of about 30 Lyman Break Galaxies 
(LBGs), most of which at 3<z<3.7, and only a few of them at 
4.3<z<5.2. A more detailed discussion on the sample selection 
is given in section 12.11 In this paper we present preliminary 
results based on a first set of data for 9 galaxies, and restricted 
to the redshift range 3<z<3.7 (Tab. [TJ. The integration times 
range from 3 to 7.5 hours on source. The goal of the SINFONI 
observations is to determine the gas metallicities by means 
of a combination of strong line diagnostics based on Hfi and 
[OIII]5007 shifted into the K band, as well as [OII]3727 and 
[NeIII]3870 shifted into the H band for sources at 3<z<3.7. 
At 4.3<z<5.2 we only rely on the [OII]/[NeIII] ratio observed 
in the K band (but sources in this redshift interval will not be 
discussed here). Petails on the gas metallicity determination are 
given in section [5] 

The two-dimensional spectroscopic capabilities of SINFONI 
are obviously exploited also to map the emission lines. However, 
the two-dimensional analysis of the spectra goes beyond the 
scope of this paper and will be discussed in a separate paper. 

2.1. Sample selection 

Galaxies in our sample are selected among z>3 galaxies pri- 
marily id entified through the Lyman-break technique, mostly 
from the Stei del et all (2003) survey and from the deep spec- 
troscopic surveys in the Chandra Peep Field South (CPFS) (e.g. 
IVanzella et aill2006l). but we also included some lensed galaxies 
(e.g. lFrve et alll2002L[2007h to better explore the low mass end 
(but none of the lensed galaxies is in the preliminary sample pre- 
sented in this paper). Galaxies were selected only amongst those 
with highly reliable spectroscopic redshift (e.g. flagged as "A" in 
IVanzella et all [2006b . We required that the redshift is such that 
the emission lines of interest for the metallicity determination 
([OIII],Hj8,[OII],[NeIII]) are out of strong sky emission lines and 
out of deep atmospheric absorption features. Actually, these re- 
quirement could not always be fulfilled for all of the emission 
lines (also because sometimes the redshift determined through 
optical spectra is not accurate, due to winds affecting UV-rest 
frame features, or IGM absorption of the Lya); however, in these 
cases any contamination by OH lines or lower S/N due to atmo- 
spheric absorption will be fully taken into account when estimat- 
ing the metallicity. 

The additional requirement is that the source has been ob- 
served with at least two of the Spitzer-IRAC bands, which at 
these redshifts sample the rest-frame near-IR light. IRAC data 
are required for a reliable determination of the stellar mass (see 
Finally, we excluded sources whose optical spectrum shows 
indications for the presence of an AGN. Moreover, for most of 
the sources we also required that deep hard X-ray data and mid- 
IR Spitzer-MIPS data are available, to better exclude the pres- 
ence of a hidden AGN, as discussed in 92. 2| 

In this paper we only present results for an initial sub-sample 
of nine sources at 3<z<3.7, for which data have been already 
obtained and reduced. The list of sources, along with some of 
their photometric properties, is given in Tab.Q] Tab.[2]lists some 
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of the physical properties of these sources as inferred from their 
broad-band spectral energy distribution. A detailed discussion 
on the extraction of these parameters will be given in $6] 

2.2. AGNs removal 

The presence of an AGN, contributing to the gas ionization, af- 
fects the observed emission line ratios. In this case the metal- 
licity diagnostic diagrams calibrated on star forming galaxies 
are not usable, since the excitation mechanism is totally differ- 
ent. As a consequence, galaxies hosting AGNs must be carefully 
avoided. 

A first step is to exclude galaxies whose optical spectrum 
(UV rest frame) shows indications for the presence of an AGN 
(e.g. NV, CIV, Hell, or broad Lya). However, the absence of 
optical-UV AGN-like lines is a required condition, but not suf- 
ficient to rule out the presence of an AGN. Indeed, even if an 
AGN is present, the associated optical-UV emission lines may 
be undetected either because obscured by dust (either on small 
scales, for the BLR, or on larger scales, for the N LR), or because 
their Narrow Line Region is not deve loped (e.g. lMaiolino et ail 
120031; iMartmez-Sansigre et aUl2006l) . 

An additional constraint comes from hard X-ray data. In 
many of the fields used by us, deep Chandra observations al- 
low the detection of obscured (Compton thin) AGNs up to z~4, 
even at Seyfert-like luminosities. Therefore, an additional re- 
quirement was that our sources are not detected in the hard X- 
rays (2-8 keV). They should not be detected also in the soft X- 
rays (0.5-2 keV) at a level higher than expected by strongly star 
forming galaxies (actually none of the galaxies is detected even 
in the soft band). In the CDFS the deep X-ray data allow us to 
exclude the presence of obscured, Compton thin AGNs with 2- 
10 keV luminosity higher than about 5 x 10 43 erg s _1 (i.e. in the 
Seyfert range). 

However, current X-ray surveys are not deep enough to 
detect Compton thick AGNs, whose emission is strongly 
suppressed even in the hard X-rays. This issue was made 
clear by recent Spitz e r resu lts (Martmez-Sansi gre et al.L [2 007; 
Alonso-Herrero et all 120061; iPolletta et all 120061; iFioreetalJ 



2007; iDaddi et all 12007). Indeed, Spitzer observations have re- 



vealed the presence of obscured AGNs, through the associated 
mid-IR hot dust emission, even in high-z galaxies that are un- 
detected in deep hard X-ray observations. These are shown to 
be high-z Compton thick (or nearly Compton thick) AGNs, 
which remained elusive to previous optical and X-ray surveys. 
Deep 24yum Spitzer-MIPS data were found to be particularly 
efficient to identify high-z obscured AGNs, even at relatively 
low intrinsic luminosities. As a consequence , we requested that 
our sources have deep MIPS data at 24/im. iFiore et ail (2007) 
showed that the mid-IR excess relative to the optical emission 
(Fi^iii/Fr), is a good tracer of obscured AGNs at high redshift. 
As listed in Tab Q] most of our sources are undetected at 24/rni. 
More specifically, all our sources have a ratio F24^ m /FR < 10, 
which is significantly lower than expected for Sy2s and QS02s 
at z~3 (F24^ m /FR ~ 20-3 10 4 ), therefore ruling out the presence 
of AGNs (both Compton thick and thin) even at Seyfert-like lu- 
minosities. 



3. Observations and data reduction 

The near-IR spectroscopic observations were obtained by means 
of SINFONI, the integral field spectrometer at VLT. SINFONI 
was used in its seeing-limited mode, with the 0.25" pixel scale 



and with the H+K grism, yielding a spectral resolution R~1500 
over the spectral range 1. 45-2.4 1/rni. 

Each target was acquired through a blind offset from a 
nearby bright star. Each observing block consisted of 10 inte- 
grations, 5 minutes each, obtained by nodding the position of 
the source within the 8" x 8" SINFONI field of view (generally 
by locating the source in two opposite corners). This observa- 
tional procedure allows background subtraction by using frames 
contiguous in time, but with the source in different locations. 
Moreover, the source was never located at the same position in 
the FOV: a minimal dithering of 0.5" was required, so that in- 
strumental artifacts can be minimized when the individual ob- 
servations are aligned and combined together. The (K-band) see- 
ing during the observations was generally better than 0.8". Each 
source was observed with a number of observing blocks ranging 
from 5 to 9. Some observing block was discarded because the 
seeing was much worse, or the background much higher, with 
respect to the other observing blocks. The total on-source inte- 
gration times are listed in Tab.Q] 

Data were reduced by using the ESO-SINFONI pipeline 
(version 3.6.1). The pipeline subtracts the sky from the tempo- 
rally contiguous frames, flat-fields the images, spectrally cal- 
ibrates each individual slice and then reconstructs the cube. 
Residual sky emission was accounted for by removing the me- 
dian of each spectral plane; this is feasible because our source 
occupy only a small part of the field of view. In some cases 
we performed an additional step in the background subtraction 
(which resulted imperfect with the previous method probably be- 
cause of minor uncertainties in the flat-fielding) by sampling the 
sky in a region outside the source (either annular or another re- 
gion in the field of view observed with the same effective inte- 
gration) and rescaling it to optimally subtract the sky lines on 
the spectrum of the source. Individual cubes where aligned in 
the spatial direction by relying on the telescope offsets and then 
averaging them together by applying a 2<x clipping to remove 
bad pixels and cosmic rays. 

The atmospheric absorption and instrumental response were 
corrected by dividing the spectrum of the scientific target by the 
spectrum of a star (spectral type OV-BV or GV) taken close to 
the source, both in time and in elevation. The intrinsic spec- 
trum of the star was removed by dividing th e observed ste llar 
spectrum by the appropriate template given in lPicklesl(fT 998). or 
by th e solar spectrum in the case of GV stars (Ma iolino et all 
1996). 



4. Results and data analysis 

We extracted the spectra within a fixed aperture of 0.75" in di- 
ameter (corresponding to ~6 kpc projected on sources at z~3.5), 
which in most cases encloses more than 70% of the emission 
line flux and generally maximizes the S/N ratio. However, one 
should keep in mind that there are metallicity gradients within 
each galaxy and therefore the aperture choice may introduce bi- 
ases, especially for what concerns the comparison with low red- 
shift surveys. This issue will be discussed more extensively in 
37.21 An exception is SSA22a-aug96M36, whose line emission 
extends significantly beyond the 0.75" aperture; in this case we 
adopted an aperture of 1.25". The resulting spectra, smoothed 
with a 2 pixel boxcar to improve the signal-to-noise, are shown 
in Figs. G][3] The location of [OII]3727, [NeIII]3870 (when ob- 
servable), H/3, [OIII]4959 and [OIII]5007 is indicated with ver- 
tical dashed lines. In each spectrum the bottom panels show the 
sky spectrum. The shaded vertical regions overlaid on each spec- 
trum highlight spectral regions affected by strong sky emission 
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CDFa-C9 
z=3.2119 




1.54 1.56 1.5S 1.60 1.62 1.64 1.66 2.02 2.04 2.06 2.08 2.10 2.12 2.14 



CDFS4414 
z=3.4714 



CDFS4417 
z=3.4733 




1.66 1.68 1.70 1.72 1.74 2.16 2.18 2.20 2.22 2.24 2.26 
*„b, [/"n] \ ab , [fj,m] 




1.66 1.68 1.70 1.72 1.74 2.16 2.18 2.20 2.22 2.24 2.26 
Kt,, [Mm] X obs [Aim] 



CDFS6664 
z=3.7967 




1.76 1.77 1.78 1.79 1.80 1.81 1.80 2.32 2.34 2.36 2.38 2.40 2.42 
[Mm] X M , [/nil] 

Fig. 1. [OIH]5007 line maps (left) and near-IR spectra (right) of the galaxies in the AMAZE sample presented here. Each [OIII]5007 
map has a size of 3" x3". The vertical dotted lines in the spectra indicate the expected location of nebular emission lines. The bottom 
panels show the sky spectrum. The shaded vertical regions overlaid on each spectrum highlight spectral region affected by strong 
sky emission lines. 



lines. In Figs.[T]|3]we also show the [OIII]5007 line map of each 
source. 

Fig.|4]shows the composite spectrum of all nine sources, ob- 
tained by shifting the spectra to the rest-frame, resampling them 
to a common wavelength scale, normalizing them by the flux of 
Hy6 and averaging them. We excluded spectral regions strongly 
affected by atmospheric absorption within individual spectra. 



The stellar continuum is detected only in a few cases, and 
even in these cases the continuum is only seen in the map pro- 
duced by stacking the cube in the spectral direction in the K or 
H band. 

The emission line fluxes were measured by fitting a sin- 
gle gaussian over a linearly interpolated, underlying continuum 
(which may be some weak stellar continuum or, more often, 
residual thermal background or residual bias subtraction). The 
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CDFS16767 
z=3.6241 




1.72 1.74 1.76 1.78 1.80 2.22 2.24 2.26 2.28 2.30 2.32 8.34 
Kb* [cm] ^„bs [cm] 



CDFS11991 
z=3.6610 





1.68 1.70 1.72 1.74 1.76 1.78 1.B02.22 2.24 2.26 2.28 2.30 2.32 
^.b, [Cm] Kb, [cm] 



CDFS2528 
z=3.6877 




1.72 1.74 1.76 1.78 2.24 2.26 2.28 2.30 2.32 2.34 2.36 2.38 
Kt,„ [cm] X oba [)im] 



SSA22a-M38 
z=3.2940 



t 




[Nell] 




1.58 1.60 1.62 1.64 1.66 1.68.06 2.08 2.10 2.12 2.14 2.16 2.: 
Kb, [cm] Kb* [cm] 



Fig. 2. Same as Fig.Q]for four additional sources. 



resulting line fluxes are given in Tab. [3] Note that for what H/3 is 
concerned we do not perform any subtraction of a stellar compo- 
nent, since the stellar continuum is always very week and gener- 
ally undetected, hence the correction for any putative stellar H/3 
is negligible. Some authors apply a fixed correction of 2A for the 
EW of a putative H/3 in absorption; in our case such a correction 
would generally affect the inferred metallicities by less than 0.03 
dex. 



5. The gas metallicity 

5. 1 . Metallicity diagnostics and calibrations 

The only method to determine the gas metallicity in faint dis- 
tant emission line galaxies is to use strong line metallicity di- 
agnostics. Essentially, the ratio between various strong, optical 
emission lines is found to depend on the gas metallicity, either 
directly and/or through other dependences (e.g. the metallic- 
ity dependence of the ionization parameter, gas density, hard- 
ness of ionizing radiation, etc.). Various strong line ratios have 
been calibrated against metallicity, either determined "directly" 
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3650 3700 3750 3800 3850 3900 4750 4800 4850 4900 4950 5000 5050 5100 
*™i [A] A„ s i [A] 

Fig. 4. Composite spectrum of the nine sources presented in this paper. 



(e.g. through the electron temperature T e method) or "indirectly" 
(e.g. through photoionization models). However, such calibra- 
tions have often been performed in relatively narrow metallic- 
ity intervals, not adequate to explore the wide metallicity range 
spanned by galaxies through the cosmic epochs, as we shall see. 
Another serious problem is that such calibrations are often in- 
consistent with each other: the same galaxy is found to have 
significantly different metallicities if different strong-line diag- 
nostics are adopted. This issue has been reviewed in detail by 
iKewlev & Ellison! d2008l) . Obviously, a wrong intercalibration 
between different metallicity diagnostics has dramatic implica- 
tions for the investigation of the metallicity evolution. Indeed, at 
different redshifts people have observed different emission lines, 
depending on the adopted band, and therefore an incorrect in- 
tercalibration between the various diagnostics may hamper the 
capability of investigating evolutionary effects, or may even in- 
troduce artificial trends. 

To minimize this issue it is recommended to use the same 
strong line calibration method for all objects at various red- 
shifts and from different surveys, or to convert the strong line 
calibrations adopted by different authors to a common cali- 
bration scale (e.g . by u sing the conversion formulas given in 
Kewlev & Ellison, 2008). However, the problem that we face in 
this paper is that no single strong line calibration method exists 
over the wide metallicity range spanned by galaxies through the 
cosmic epochs (as we shall see). Some methods nominally span 
a somewhat wider metallicity range, but they are known to run 
into serious troubles in some metallicity intervals. 

The electron temperature T e method (e.g. by exploiting 
the intensity of the [OIII]4636 auroral line) provides a good 
measure of the metallicity below about 12+log(0/H)<8.3 (e.g. 



iPilvuginl 1200 U iPettini & Pagel |2004|) . and can be used to cal- 
ibrate the strong line ratios in this range. The reliability of 
the T e method in the low metallicity range is confirmed by 
the comparison with the stell a r (OB) photos pheric metallicity 
measurements dBresolin et all l2006bl l2007bl) . The T e method 
has been extended to higher metallicities by various authors 
dKennicutt et all 120031: iGarnett et all 120041: iLiang et all 120071: 
lYin et all 120071) . However, at high metallicities the T e method 
tends to saturate and to underestimate significantly the true 
metallicity, due to temperature fluctuations and gradients, both 
within individual HII regions a nd over the whol e galaxy. This 
issue is expected theoretically ( Stasiriska, 2005) and verified 
observationally by the comparison with the metallicities deter- 
mined through recom bination l i nes, which are insensitive to tem- 
perature fluctuations (Bresolin, 2006a, 20071). 

Photoionization mo dels are an alternative way of calibrating 
strong line ratios (e.g . iTremonti et aTl l2004t IKewlev & Dopita, 
2002; Zaritsk vet all .11994 . especially at high metallicities, 
where most of these studies apply. However, all photoionization 
models are subject to significant uncertainties and possible sys- 
tematic effects. The observed spread i n calibration between dif - 
ferent models highlights this problem dKewlev & Ellison. 2008). 
The photoionization models presented in IKewlev & Dopita 
d2002l) are probably not free from the uncertainties and possi- 
ble systematic effects discussed above, however they provide re- 
sults which are intermediate among other photoionization mod- 
els dKewlev & Ellisonl [2008), and therefore can be considered 
fairly representative of this class of calibrations. Moreover, in- 
dependent "direct" determinations of the metallicity (by exploit- 
ing the temperature-ins ensitive method of recombination lines, 
lBresoUdl2006aLl2007ah are in fair agreement with the photoion- 
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Fig. 5. Rela t ions b etween strong emission line ratios and gas metallicity. Blue squares are low metallicity galaxies (from 
Nag ao etail 120061) for which the metallicity is inferred through the electron temper ature T e method. Black dots are SDSS DR4 
galaxies for which metallicities are inferred through the photoionization models in iKewlev & Dopital ((2002). The red line curve 
shows a polynomial fit to the relations in the form given by Eq.Q]with the coefficients given in Tab. [8] 



ization models provided by IKewlev & Dopital ([2002). The lat- 
ter did not investigate photoionization models at metallicities 
12+log(0/H)<8.4. Other studies attempt to extend photoioniza- 
tion models to 12+log( Q/H)<8.3, but fai l to reproduce the ob- 
served line ratios (e.g. lDopita et al.Ll2006l) . 

Since no single method is capable of providing a calibration 
of the strong line diagnostics over the wide metallicity range re- 
quired to sample the evolution of galaxies through the cosmic 



epochs (7.7<12+log(0/H)<9.1), in this paper we have to com- 
bine two different methods depending on the metallicity range. 

At low metallicities (12 + log(0/H) < 8.35) we use cali- 
brations of the strong-line diagnostics based on the T e method, 
which is regarded as reliable at low metallicities and free 
of problems rela t ed to temperat u re fluctuations and gradients 
dStasinskal 120051: iBresolin et all l2006bl l2007bl) . We use the 
sampl e of 259 low metallicities galaxies gathered bv lNagao et alj 
(2006) for which measurement of the auroral line [OIII]4636 is 
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Fig. 6. Example of the diagnostic tools used to determine the metallicity in the specific case of the composite spectrum. The upper 
left panel shows the best solution (blue cross) and the Icr confidence level in the Ay-metallicity plane. In the other panels the black 
solid line (best fit) and the dashed lines (dispersion) show the empirical relations between various line ratios and the gas metallicity 
(see Fig . O . The green errorbars show the observed ratios (along the Y-axis) and the best-fit metallicity with uncertainty (along the 
X-axis); the blue cross shows the de-reddened ratios, by adopting the best-fit extinction; the red line shows the projection of the Icr 
uncertainty of the fit in the top-left panel. 



available. Th e metallicity determ ination based on the T e method 
is detailed in Naga o etail d2006l) . The inferred metallicities for 
this sample are in the range 7.1 < 12 + log(0/H) < 8.3. 
The inferred empirical relations between metallicity and various 
strong-line diagnostics are shown in Fig. [5] with blue squares. 
Note that the T e method is only used to calibrate the strong line 
diagnostics in local galaxies (and at low metallicities), but it will 
not be used to directly measure metallicities in our high-z ob- 
jects, since the [OIII]4636 line is too faint to be detected. 

At 12 + log(0/H) > 8.35, where the T e method is know 
to fail (as discussed above), we have to rely on photoionization 
models. Mor e specifically, we adopt th e calibrations provided by 
the models in Kewley & Dopita (2002), aware of the caveats dis- 
cussed above, and in particular that these models are valid only 
at 12 + log(0/H) > 8.4. We then used the data of star form- 
ing galaxies in SDSS DR4, by adopting the s ame c onstraints on 
the signal-to-noise discussed in lNagao et alj d2006t) . and by also 
considering only o bjects with log( r NII 1 6564 /rOII137271) > 1.2, 
as recommended bv lKewlev & Dopita (2002). This selection re- 
sults in a total of 22,482 objects. After determining t he gas 
metallicity for each object with the iKewlev & Dopital (2002) 
method, we derived the empirical relations with various strong- 
line diagnostics as shown in Fig.[5](black dots). 

The relations obtained by combining both the low metallicity 
and the high metallicity samples were fitted with a polynomial 
curve. To avoid the fit to be dominated by the regions contain- 
ing the largest number of objects (i.e. by the SDSS sample at 
8.7 < 12 + log(0/H) < 9.1), we divided the relations in metal- 
licity bins (generally spaced by 0.1 dex), we derived the me- 
dian and dispersion of the strong-line flux ratio and of the stellar 
masses within each bin, and then fitted the polynomial function 



to these medians. In most cases a second order or a third order 
polynomial is appropriate to describe the relation over the full 
metallicity interval. However, in some cases a fourth order poly- 
nomial is required. The general functional form for describing 
the strong-line metallicity calibration is therefore: 

log/? = Co + c\x + C2X 2 + C3JC 3 + C4X 4 (1) 

where log R is the logarithm of the strong-line ratio, and x is the 
metal licity relative to sola r (x = log (Z/Z ) = 12 + \og(0/H) - 
8.69, lAllende Prieto et al.1 d2001l) ). The coefficients c - c 4 for 
each strong-line ratio are listed in Tab. [8] and the resulting best 
fit polynomials are plotted with a solid red curve in Fig. [5] 

5.2. Determination of the gas metallicity 

In the redshift range 3<z<3.7 investigated in this paper 
[OII]3727 and (not always) [NeIII]3870 are observable in the H 
band, while H/3 and [OIII]5007 are observable in the K band. 
This allows us to use five metallicity diagnostics (not all of 
them independent of each other), namely: R23 = ([OII]3727 + 
[OIII]4959 + [OIII]5007)/Hje, [OIII]5007/Hy8, [Oil] 3727 /H/3, 
[OIII]5007/[OII]3727, [NeIII]3870/[OII]3727. The relation- 
ship between these ratios and the gas metallicity, along with 
their dispersion, are shown in Figs. [5] and [6] Each of these 
diagnostics has advantages and disadvantages. For instance, 
[OIII]5007/Hj6 is essentially unaffected by dust reddening, but 
it has a double metallicity solution for each value of this ratio. 
[OIII]5007/[OII]3727 has a monotonic dependence on metal- 
licity, but it is potentially affected by dust reddening and it is 
also affected by a larger dispersion. [NeIII]3870/[OII]3727 is 
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both a monotonic function of metallicity and little affected by 
dust reddening, nonetheless [Nelll] is generally the faintest of 
all these lines and the most difficult to detect. However, if these 
various diagnostics are used simultaneously, then it is possible 
to both account for dust reddening and remove the ambiguity of 
double solutions. Essentially, if all the diagnostics are used then 
only some combinations of metallicity and dust reddening are 
allowed by the data. 

More specifically, we selected the following independent 
metallicity diagnostics: [Offl]5007/Hj6, [OIII]5007/[OII]3727 
and (when available) [NeIII]3870/[OII]3727 (note that the ratio 
[OIII]5007/Hy6 is very similar to the "classical" R23 parameter 
at the low metallicities investigated by us, since at low metallic- 
ities [OIII]/[OII] is high). We then determined the best pair of 
metallicity and extinction that minimizes the x 1 in the three cor- 
responding diagrams, both by including the measurement errors 
and the dispersion of each calibration diagram. In practice, since 
[Nelll] is often undetected or not observable, the metallicity is 
mostly determined through the [OIII]5007/HjS ratio, while the 
[OIII]5007/[OII]3727 ratio is used to discriminate which of the 
two metallicity solutions for the [OIII]5007/Hy6 ratio applies and 
also to provide some constraints on the dust extinction (although 
the latter has generally negligible effect on the metallicity deter- 
mination since [OIII]5007/Hy8 is insensitive to reddening). Q 

As an example, we show the results of this method in Fig. [6] 
in the case of the composite spectrum. The upper left panel 
shows the best solution (blue cross) and the lcr confidence level 
(red curve, obtained from solutions with Ax 2 = 1) in the Ay- 
metallicity plane. In the other panels the black solid line (best 
fit) and the dashed lines (dispersion) show the empirical relations 
between various line ratios and the gas metallicity (Fig. |5); the 
green errorbars show the observed ratios (along the Y-axis) and 
the best-fit metallicity with uncertainty (along the X-axis); the 
blue cross shows the de-reddened ratios, by adopting the best-fit 
extinction; the red line shows the projection of the lcr uncer- 
tainty of the fit obtained in the top-left panel. It can be noted that 
the extinction is subject to a large uncertainty, but the metallic- 
ity is relatively well constrained. The metallicities resulting from 
the procedure discussed above are listed for all objects and for 
the composite spectrum in Tab. [3] 

5.3. Evolution of the metallicity diagnostics 

One of the main worries when using strong emission line diag- 
nostics in high-z sources is that the empirical calibrations are ob- 
tained by using local sources. Since the dependence of the strong 
lines ratios on metallicity involves also other dependences (e.g. 
on the ionization parameter, on the shape of the ionizing radi- 
ation, on the gas density), the evolution of the average galaxy 
properties (e.g. SFR, compactness) with redshift may affect the 
calibration of the metallicity diagnostics. It is very difficult to in- 
vestigate this issue, since in principle one would need to obtain 
an empirical calibration at high-z by observing primary metallic- 
ity tracers (e.g. [OIII]4636), which are however extremely faint. 
An alternative way is to construct diagrams which are sensitive 
to the excitation mechanism, by disentangling the dependence 
on metallicity, and verify whether the line excitation conditions 
change with redshift. 

Such a test was performed, at lower redshifts, by investigat- 
ing the diagram [OIIIJ/H/? versus [NII]/Ha for sources where 



1 For the Ay -meta llicity fit we used the Milky Way extinction curve 
dCardelli et all IT9 89j) (Rv = 3.1), which is often preferred for the neb- 
ular lines. 



all of these lines are ob servable dShaplev et all 120051 ; lErb et ail 
20061 iLiu et all 120081) . It is found that a fraction of sources at 
z~l-2 are offset with respect to the sequence described by local 
HII galaxies, and displaced towards the AGN locus. An interpre- 
tation is that at leas t some of t hese sources are affected by some 
AGN contribution (ILiu et al ., 2008), which were not excluded 
from the sample because elusive in the UV rest-frame spectra. 
Other sources may be characterized by truly different physical 
conditions with respect to local HII regions, and in particular 
high er ionization parameter. How ever, even in the latter cases 
both lBrinchmann et all (2008) and lLiu et ail d2008l) find that the 
calibration of the strong line metallicity diagnostics do not de- 
viate by a large amount with respect to local HII galaxies. In 
particular, they find deviations by only about 0. 1 dex (or less) in 
terms of metallicity calibration, depending on the specific diag- 
nostic adopted. 

For what concerns the sources at z~3.5 in AMAZE, the 
available emis sion lines allow us to construct the so-called 
BPT diagram dBaldwin et all fl98~ll) . i. e. rOIII15007/H8 v ersus 
[OIII]5007/[OII]3727. As discussed in iDopita et all (120061) this 
diagram is strongly degenerate in terms of metallicity, but it is 
sensitive to both the ionization parameter and the hardness of 
the ionizing source. Within the observational uncertainties, the 
sources in our AMAZE sample do not deviate from the sequence 
of local HII galaxies on the BPT diagram, suggesting that the 
excitation conditions do not differ significantly from the local 
galaxies used to calibrate the strong metallicity diagnostics. 

6. Stellar masses 

To derive the stellar masses for the LBGs in the AMAZE sample, 
we used an approach based on broad-band spectral fitting tech- 
nique. Broad-band photometric data for the sources in the CDFS 
were collected from the G OODS-MUSIC multiwavelength cat- 
alog (Gra zian et all 120061) . This catalog provides photometric 
data in 14 spectral bands (from UV to the Spitzer-IRAC bands), 
and it has been recently upda ted to include the S pitzer-MIPS 
data at 24/zm. For the LBGs in lSteidel etai] d2003l) . optical pho- 
tometric data ( U,G,R,I) were ex tracted from the publicly avail- 
able images dSteidel et alll2003l) . while Spitzer IRAC and MIPS 
data were obtained from the Spitzer archive; the photometry ex- 
tr action was performed following the same methods described 
in lGrazian etaT] d2006l) . 

The SED fitting technique adopted here is the same as in pre - 
vious papers (Fontana et al., 2006; Grazia n et al 1 120061 12007). 
and s imilar to those adopted by other g r oups in the literature 
(e.g iDickinson et all 120031; iDrorv et all l2004t iPozzetti et all 
2007). This technique is based on comparing the observed mul- 
ticolor distribution of each object and a set of templates, com- 
puted with standard spectral synthesis models (see below) and 
chosen to broadly encompass the variety of star-formation his- 
tories, ages, metallicities, and extinction of real galaxies. More 
specifically, we considered exponentially decaying SFR with e- 
folding times ranging from 0.1 to 15 Gyr. We used the Salpeter 
IMF (M„„„ = 0.1 M and M mclx = 65M Q ), ranging over a set of 
metallicities (from Z = O.O2Z to Z - 2.5Z ) a n d dus t extinc- 
tion (0 < E(B -V) < 1.1, with a lCalzettietail d2000l) attenu- 
ation curve, which is generally more appropriate for the stellar 
component). For each model of this grid, we computed the ex- 
pected magnitudes in our filters set and found the best-fitting 
template with a standard^ 2 normalization. The stellar mass and 
other best-fit parameters of the galaxy, like SFR, age, and dust 
extinction, are fitted simultaneously to the actual SED of the ob- 
served galaxy. The metallicity of each galaxy is fixed to value 
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Fig. 7. Mass-metallicity relation observed at different r edshifts. The blue, solid a nd dotted lines indicate the mass-metallicity relation 
and its dispersion observed at z~0.07, as inferred by iKewlev & Ellisonl d2008l) . The red diamonds with errorbars show the mass- 
metallicity relation at different redshifts traced by individual objects (or by stacked spectra in the case of panel b, z~2.2). In panel c 
(z~3.5) the green square with dashed errorbars is the composite spectrum (which has been assigned the median mass of the sample, 
but slightly offset to avoid a confusing overlap with an individual source). Black, solid lines show the analytical function Eq.[2]with 
the best fitting parameters listed in Tab. [5] 



closest to the one determined by us through the nebular lines 
(Tab.|3). 

The stellar mass derived here is subject to uncertainties and 
biases related to the synthetic libraries used to carry out the fit- 
ting of the galaxy SEDs. In general, the stellar mass turns out 
to be the least sensitive parameter to variations of the input 
model assumptions, and the extension of the SEDs to mid-IR 
wavelengths (near-IR rest-frame) with IRAC greatly reduces the 
formal uncert a inties on the derived stellar masses, as shown in 
Font ana et al.l d2006l) . The uncertainties in the stellar mass are 
derived as follows: we compute the 90% confidence level on the 
mass by scanning the^ 2 levels, fixing the redshift and the metal- 
licity for each galaxy but allowing the other parameters (SFR, 
age, dust extinction) to change. 

Age and star formation rate are more uncertain parameters 
to derive. In some cases we formally obtain best-fit ages be- 
low ~50 Myr, which are below the dyn amical timescales fo r 
the star forming regions in these systems (Shap lev et all 12001). 
Moreover, the conversion between UV luminosity and SFR be- 
comes highly non-linear below this age. As a consequence, we 
decided to restrict the allowed ages to >50 Myr. However, this 
choice may only affect the inferred SFR, while the determination 
of the stellar mass is essentially unaffected, as discussed above. 

For what concerns the library o f spectral synthesi s mod- 
els we adopt both those provided by Bruzual & Chariot] (2003) 
(hereafter BC03) and those by Marastonf ?2005l) (hereafter M05). 
The resulting stellar masses are tabulated for both cases in 
Tab. |2] The masses inferred by using the M05 models are pre- 
ferred, since they take into account the contribution by TP- 
AGB stars, and may differ from the stellar masses obtained with 
BC03 by even a factor of two, especially in older stellar sys- 
tems. However, previous works on the mass-metallicity relation 
at lower redshift adopted the BC03 templates. Therefore, when 
comparing our results with the previous works at lower redshifts 
we will adopt for consistency the masses obtained with the BC03 
templates. In Tab. [2] we also list the SFR, age and reddening in- 
ferred by the SED fitting (adopting templates by BC03, again for 
a consistent comparison with previous works). 



7. The mass-metallicity relation at high redshift 

7.1. Comparing the mass-metallicity relation at different 
redshifts 

Different studies of the mass-metallicity relation at various 
redshifts have employed different diagnostic lines and differ- 
ent calibrations. As disc ussed in and more extensively in 
IKewlev & Ellisonl (2008), the mismatch between the different 
calibration scales may introduce artificial evolutionary effects of 
the mass-metallicity relation. Therefore, it is important that dif- 
ferent strong-line diagnostics used in different surveys are cross- 
calibrated in a consistent way. The relations obtained in J5]p r< > 
vide such a common cross-calibration between different strong- 
line diagnostics on the same metallicity scale. In this section we 
apply (when required) the correction to the metallicities inferred 
by past surveys at lower redshift to match our metallicity scale. 
We also apply corrections to the mass scale to account for the 
different IMF's adopted by previous works. 

As discussed in jjQ] the local (z~0.1) mass-metallicity rela- 
tion was derive d bvlTremonti et al.l d2004l) by using SDSS spec- 
tra from DR2. IKewlev & Ellisonl d2008l) re-determined the lo- 
cal mass-metallicity relation by using SDSS spectra from DR4 
by setting tighter limits on the redshift range (0.04<z<0.1) so 
that the projected SDSS fiber covering factor is >20% of the 
total photometric g'-band light, and also to minimize incom- 
pleteness effects at higher redshi fts. The resulting media n red- 
shift of their sample is -0.0 7. IKewlev & Ellison! d2 008) cali- 
brate the metallicities with the IKewlev & Dopital d2002h method, 
which is the same adopted by us at 12+log(0/H)>8.35, hence 
no additional correction is required to match our metallicity 
scale . The only correction to apply is for the stellar m asses, 
since lTremonti et all d2004b andlKewlev & Ellisonl d2008l) adopt 
a different IMF dKroupal l200lh . We calculate that the masses 
in IKewlev & Ellisonl d2008l) must be multiplied by a factor of 
1.17 to comply with our IMF (note that the IMF's differs not 
only in terms of shape but also in ter ms integration limits). The 
thin blue solid lines in Fig.|7]show the IKewlev & Ellisonl (2008) 
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mass-metallicity relation corrected as discussed above. The blue, 
dotted lines indicates the lcr dispersion of the same rela ti on. 

At 0.4<z<l we use the results by ISavaglio et al.l d2005l) . 
For consistency with our calibration scale we re-determine the 
metallicities for each object in their sample by applying the 
same procedure described in 35.21 to the line fluxes tabulated 
by them. We exclude from their sample objects without K-band 
data, since in these cases the stellar mass uncer tainties are too 
large. We also have corrected the stellar masses in Savagli o" et alj 
(2005) by a factor of 1.4 to comply with the IMF adopted by us. 
The resulting mass-metallicity relation at z~0.7 is shown with 
red diamonds and errorbars in Fig. [TJi. 

At z~2.2 we use the results by lErb et al.l d2006l) . who infer 
the metallicity of LBG's at this redshift through the [NII]/Ha 
ratio measured in stacked spectra. Also in this case we re- 
determine the metallicity in each mass bin by using the [NII]/Ha 
metallicity calibration obtained in ^5] to be consis tent with 
the ca librations adopted by us. The stellar masses in Er b et alj 
(2006J) have to be corrected by a factor of 1.4 to comply with 
the IMF adopted by us. The resulting mass metallicity relation 
at z~2.2 is shown with red diamonds and errorbars in Fig.|7}5. 

Finally, Fig. [7]: shows the mass-metallicity relation inferred 
from the initial sample of nine AMAZE sources at z~3.5. The 
red diamonds with solid errorbars are individual objects. The 
green square with dashed errorbars is the composite spectrum. In 
this plot, for consistency with the other works at lower redshifts, 
we use the stellar masses inferred with the BC03 templates. 

Previous versions of our mass-metallicit y relation z~3 , 
presen ted in our previous preliminary works ( Mai olino et al.L 
l2007allbl) . were slightly different because of lower S/N spectra 
and also because we used different calibrations (both for metal- 
licity and stellar masses). 

For a more straightforward comparison of the mass- 
metallicity relation at different redshifts, it is useful to describe 
these relations by fitting them with the same functional form. In 
order to minimize the number of free par ameters we find statis - 
tically satisfactory an approach similar to Sava glio"et al.l (|2005): 
the quadratic function fitting the local mass-metallicity relation 
is shifted in mass and in metallicity to provide the best fit of the 
mass-metallicity relation at various redshifts. More specifically 
we adopt the following description of the mass-metallicity rela- 
tion: 

12 + log(0/H) = -0.0864 (logM, - logM ) 2 + K (2) 

where logMo and Ko are determined at each redshift to obtain 
the best fit to the observed data pointfl as listed in Tab. [5] 

7.2. Aperture effects 

Since galaxies, and especially disk galaxies, are often character- 
ized by metallicity gradients (the metallicity decreasing towards 
the outer regions), a possible caveat when comparing metallici- 
ties at different redshifts is the different aperture projected on the 
source. In particular, at high redshift spectroscopic observations 
are likely to sample most of the galaxy, while at low redshift 
and in local galaxies the spectroscopic aperture samples mostly 
the central higher metallicity region. This effect may mimic a 
metallicity evolution. 

When comparing surveys at high redshifts this should be a 
minor issue, since the projected apertures on the sources are not 

2 Note that, to avoid the fit to be dominated by the mass intervals with 
largest number of objects, we first obtained metallicity medians within 
mass bins and then fitted the resulting values with Eq.[2] 



very different. At z~ 0.7, the adopted aperture of 0.75-1.3" in 
Sava glio et al.l (120051) corresponds to about 5.6-9.3 kpc; at z~2.2 
lErb et al.l d2006l) adopt an aperture of 0.76" corresponding to 
about 7.2 kpc. Our aperture of 0.75" at z~3.5 corresponds to 
about 6 kpc. 

However, aperture effects may be more serious for the lo- 
cal sample. At the median redshift of 0.07 the SDSS fiber size 
(3") has a median size of 4 kpc, and a median covering factor 
~34% relative to the total g'-band light, in contrast with a cov- 
ering factor of ~70% at z~3.5. The aperture effect is stronger 
for local high mass galaxies, which are generally bigger and 
for which the covering factor reach values as low as 20%. The 
problem of a differential aperture effect as a function of galaxy 
mass may also affect the shape of the local mass-metallicity re- 
lation, making the obser ved relation steeper than it actually is 
dKewlev & Ellisonl 120081) . However, the absolute magnitude of 
this effect in the local sample is e stimated to be at most ~0.1- 
0.15 dex dKewlev & Ellisonl [2008). which is significantly lower 
than the metallicity evolution observed at high redshift, at least 
at z~3.5. 

Additional issues related to aperture effects will be discussed 
in the context of the comparison with models in 37.51 

7.3. Selection effects 

When comparing the mass-metallicity relation at z~3.5 inferred 
from LBG's with that inferred from local (or lower redshift) 
samples of star forming galaxies, one must be aware that we are 
comparing different classes of objects, which are not necessarily 
linked from an evolutionary point of view. As a consequence, the 
evolution of the mass-metallicity relation inferred in this paper 
should be regarded as the evolution of the mass-metallicity re- 
lation of galaxies representative of (or contributing significantly 
to) the density of star formation at each epoch, and not the evolu- 
tionary pattern of individual galaxies. This issue will be further 
discussed in the next sections. In this section we mostly inves- 
tigate whether galaxies in our sample are representative of star 
forming galaxies at z~3.5. 

Stellar masses of LBGs at z~3 were measured by 
Shap lev et al.l d200llP l but for a subsample of galaxies about 
one magnitude brighter (in R-band) than the parent sample of 
ISteidel et ail d2003l) . IShaplev et ail (1200 ll) obtain a median stel- 
lar mass of 2.4 x 10 10 M . The galaxies in our sample (Tab.|2]and 

have a median stell ar mass of ~ 1.6 xl0 1 M Q , which is close 
to the value obtained by Shap lev et al.l d200 1 ) for their large LBG 
sample. In any case, the distribution of stellar masses is not a 
concern, since the mass is one of the two variables that we are 
mapping on the mass-metallicity relation: even if we had a bias 
in terms of stellar mass, this would simply imply that we pref- 
erentially populate the diagram in a certain mass range, making 
the estimation of the mass-metallicity relation more uncertain in 
other mass ranges (because under-populated), but not biased. 

Biases in terms of star formation rate are of a greater con- 
cern. Our selection requirement that sources must have a highly 
reliable spectroscopic redshift may bias our sample towards 
sources with strong UV continuum or strong Lya, hence higher 
than average SFR. Tab. [2] shows the SFR inferred from the rest- 
frame UV continuum of the sources in our sample (see $6] for 
details), from which we infer a median SFR of ~ 100 M Q yr _1 . 
This is similar to the median SFR (90 M yr~') of LBGs at 

3 Note that lShaplev et alj d200lh did not have access to Spitzer-IRAC 
data, and therefore their masses may be subject to significant uncertain- 
ties. 
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z~3 obtained by Shap lev et alj d2001l) (who adopted a simi- 
lar approach as ours to estimate the SFR from the UV contin- 
uum). However, since the latter work is biased towards slightly 
brighter optical magnitu des, the median SFR of the LBGs at z~3 
in the parent sample of Steidel et al. (2003) is probably some- 
what lower. Since the median stellar masses are similar, then 
the specific star formation rate (SSFR) in our sample may be 
some what higher (up to a factor of 2) than in the LBG sam- 
ple of lSteidel etail d2003l) . lEUison et all d2008l) investigated the 
effect of the SSFR on the metallicity in local galaxies. They 
found that, for a variation of the SSFR by a factor of two, the 
metallicity varies by less than ~0.1 dex in low mass galaxies 
(M, < 1O 1O M ), while no metallicity variations are found in 
massive galaxies (M* > 1O 1O M ). As a consequence, a possible 
bi as of our sample in terms of SSFR relative to the LBG sample 
of lSteidel et alj d2003l) should not affect significantly the inferred 
mass-metallicity relation. 

Another possible source of bias is that LBGs are selected 
through their UV rest-frame colors, hence they mis s any pop- 
ulatio n of heavily reddened star forming galaxies. iReddv et alj 
d2007l) estimate that, due to this selection effect, LBG's rep- 
resent ~ 47% of the population of star forming gal a xies a t 
z~3 with R < 25.5 (see also iHopkins & Beacoml [2006). 
Dust reddened galaxies are naively expected to be more metal 
rich ( since metallicity a nd dust content cor relate, iHunt et all 
120051) . On the contrary, iRupke et all d2008l) have shown that, 
at least at low-z, dusty galaxies (IR-selected) are characterized 
by gas metallicities lower than optical- and UV-selected galax- 
ies (probably due to metal poor gas infallin g from the outskirts 
in merging/interacting systems). Recently Cap uti et"ai] ([2008) 
have found a similar effect in dusty galaxies at intermediate red- 
shift (i.e. metallicities lower than in optically selected galaxies). 
High-z dusty objects may behave differently. However, if the 
same phenomenon is also present at high-z, then LBG's would 
provide an upper limit to the metallicity of galaxies at z~3. 
However, we do not speculate further on the properties of star 
forming galaxies that are not sampled by LBGs, and we sim- 
ply emphasize that the results presented in this paper only apply 
to about half of the star forming galaxies at z~3, i.e. those UV- 
selected. 

7.4. The evolution of the mass-metallicity relation 

The evolution of the mass-metallicity relation is summarized in 
Fig. [8] where we plot the best fits resulting from Eq. [2] both at 
z~3.5 from the AMAZE survey, and at lower redshifts from pre- 
vious surveys. 

Figs-LZHHJhighlights a clear evolution of the mass-metallicity 
relation of star forming galaxies through the cosmic epochs. As 
already discussed in jj7.3l this evolution should not be seen as 
the evolutionary sequence of individual objects, since at each 
redshift the various surveys are sampling different classes of star 
forming galaxies, which are not necessarily each other progeni- 
tors. The trend observed in Figs. [7|-[8] should be regarded as the 
evolution of the mass-metallicity relation of galaxies dominating 
(or contributing significantly) the star formation density at each 
epoch. 

At M, ~ 1O IO M the metallicity at z~2.2 is lower by a fac- 
tor of about 2.5 with respect to local galaxies. Even if highly 
significant, such metallicity decrease is modest if one considers 
that from z=0 to z=2.2 the elapsed time is ~11 Gyr, i.e. about 
75% of the age of the universe. From z~2.2 to z~3.5 the aver- 
age metallicity of galaxies decreases by another factor of about 

2.5. However, the latter evolution is much stronger, and faster. 
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Fig. 8. Comparison of the mass-metallicity relation observed at 
different redshifts, as parametrized by the analytical function 
Eq.|2]and coefficients in Tab. [5] 
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Fig. 9. Average metallicity of star forming galaxies as a function 
of the cosmic age of the Universe, relative to local galaxies, for 
three different families of galaxies with different stellar masses 
(M»/M = 10 9 , 10 10 , 10 11 ). Ages are sampled at the redshifts 
where observations are directly available. The average metallic- 
ity for each stellar mass was inferred by using the parametriza- 
tion given by Eq.|2]and Tab. [5] 



Indeed, such a metallicity variation occurs on a much shorter 
time scale, only ~1 Gyr. This effect is shown in more clearly 
Fig-E| where the the average metallicity of star forming galaxies 
is plotted as a function of the age of the universe, for different 
classes of galaxies with different stellar masses, by exploiting 
the analytical function in Eq. [2] (note that, as discussed above, 
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Fig. 10. Comparison between models/simulations predictions for the mass-metallicity relation at z~3 and the results of our obser- 
vations. The predictions obtained by models and simulations are shown with green points and dashed lines (the aperture used in 
the models/simulations is reported when available). The solid, black lines show the mass-metallicity relation observed at the same 
epochs (interpolated to the same, exact redshifts of the simulation by using Eq.[2]and Tab. [5]). The thin black line shows the observed 
mass-metallicity relation by using the masses estimated with the BC03 templates, while the thick line is the inferred relation when 
using the M05 templates. 



this figure does not provide the evolution of individual galax- 
ies). We further note that the evolution is strong even in massive 
galaxies. Clearly, z~3.5 is an epoch of major action for the evo- 
lution of galaxies, both in terms of star formation and chemical 
enrichment, even for massive systems. 

We note that a simil ar strong evolution o f the metallicities 
at z>3 was obtained by iMehlert et al.1 d2006l) by investigating 
the stellar metallicities of a few galaxies with bright UV contin- 
uum. Their result was however affected by uncertainties on the 
absolute calibration of the stellar metallicity tracers. It is also 
interesting to note that the metallicities obtained by us at z~3.5 
are in fair agreement with those expected at the same redshift 
by iPanter et ail d2008l) . who inferred the metallicity evolution 
of galaxies by modelling its "fossil" spectral signatures in local 
galaxies. 

The additional interesting result is the indication of a differ- 
ential, mass-dependent evolution of the metallicity. In particular, 
the metallicity evolution in low mass systems appears stronger 
than in massive galaxies (Figs. [7]-[9j>. This finding requires more 
statistics to be confirmed at z~3.5 (to come with the completion 
of the AMAZE project). However, the evolution of the slope of 
the mass-metallicity relation (relative to the local slope) is signif- 
icant even at z~0 . 7 (as already noted by ISavaglio et all 12005; 
Kobu lnickv et all 12003) and at z~2.2. A detailed investigation 
of differential selection effects as a function of stellar mass is re- 



quired to rule out that observational biases are not affecting the 
slope of the mass-metallicity relation at each epoch. However, if 
confirmed, such mass-dependent evolution of the metallicity can 
be regarded as the "chemical" version of the galaxy downsizing: 
high mass galaxies reach high metallicities at high redshift, on 
short timescales, while low mass systems enrich their ISM over 
a prolonged period of time, extending to the current epoch. 

7.5. Comparison with models of galaxy evolution 

As mentioned in the introduction there is an intense theoretical 
activity aimed at interpreting the nature and origin of the mass- 
metallicity relation, and also at providing predictions on the ex- 
pected mass-metallicity relation at high redshift. Comparing the 
models predictions with the observational results in a consistent 
way is not simple. Indeed, theoretical models predict a variety of 
galaxy populations, spanning a wide range of properties (e.g. in 
terms of SFR), while observations are limited to samples match- 
ing the survey selection criteria. Moreover, generally theoretical 
works provide the metallicity integrated over large apertures, in- 
cluding most of the galaxy, while observational metallicity mea- 
surements are generally obtained within a smaller aperture. In 
the future a collaborative effort with various theoretical groups is 
planned to match the outcome of models and simulations to the 
observational selection effects. However, a preliminary compar- 
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ison with the already published theoretical results is instructive 
to in fer some initial cons traints on galaxy evolutionary models. 

Savaglio et al. (2005) interpreted the evolution of the mass- 
metallicity relation from z~2.2 to z=0 through a closed-box 
model with an exponentially decaying SFR oc e~'/ T , where the e- 
folding time r decreases as a function of galaxy mass. However 
the inclusion of our data at z~3.5 makes this model not suit- 
able. It is difficult to simultaneously fit the observed mass- 
metallicity relations at z=0.1, z=0.7, z=2.2 and z=3.5 with a 
simple closed-box model, unless more complex scenarios of the 
SF history are envisaged. Moreover, an e xponentially decaying 
SFR, with the e-folding times provided by Savaglio et al. (2005) 
makes the SFR extrapolated to local massive systems well below 
1 M yr , i.e. these should be local quiescent galaxies (probably 
massive elliptical), which cannot be representative of the local 
star formi ng galaxies used to de rive the local mass-me t allicit y 
relation in lTremonti et all d2004l) and lKewlev & Ellison! (120081) . 
Finally, it is unlikely that the closed-box scenario applies to 
LBGs, which are c haracterized by strong, unbound winds (e.g. 
iPettini et all 120021) . 

Withi n the framework of th e hierarchical models of galaxy 
evolution, Ide Rossi et al.l d2007l) performed numerical hydrody- 
namical simulations enabling them to provide detailed predic- 
tions on the evolution of the mass-metallicity relation at vari- 
ous epochs. In the top-left panel of Fig. [10] the dashed green 
line shows the mass-metallici t y rel ation predicted at z=3 by the 
simulations of Ide Rossi et al.l d2007l) (therein Tab. 3), while the 
shaded areas give t he dispersion i n ferred from the same sim- 
ulations (note that de Ros si et all 120071 use the same IMF 
adopted by us, therefore no further correction is required). Since 
we do not have the observational data at exactly the same epoch 
of the simulations (z=3), we interpolate the observed mass- 
metallicity to the same epoch of the simulations by using Eq. [2] 
and Tab. [5] The observed mass-metallicity relation at z=3 is 
shown in Fig. [10] with a thin, black solid line. We also show the 
mass-metallicity relation inferred by adopting the stellar masses 
measured with the M05 templates (thick solid line). Fig. [10] 
shows a significant discrepancy between simulations and obser- 
vations. The discrepancy was also presen t at z- 2 (although at 
a lower level), as noted by Ide Rossi et alj (|2007). They suggest 
that the inconsistency between simulations and observations is 
due to the lack of significant SN feedback in their simulations, 
which would remove metal enriched gas and lower the global 
gas metallicity. 

However, the discrepancy with observations is present also 
for simulations that include the effect of SN feedback. This is 
the case of the hierarchi cal, three-dimens i onal c hemodynamical 
simulations presented in Kobavas hi et al.l (|2007), which include 
the feedback from SNII and hypernovae. They provide predic- 
tions of the mass-metallicity relation as a function of redshift. 
The green points in the top-right panel of Fig. [10| show the re- 
sults of such simulations at z=3 (also Kobavashi et al., 2007, 
use the same IMF adopted by us, therefore no further correc- 
tion is required). Even in this case there is a discrepancy be- 
tween simulations and obser vations. The discrep ancy is much 
reduced with respect to the Ide Rossi et alj d2007l) simulations, 
but still significant at low masses (log(M*/M ) < 10). At high 
masses (log(M»/M ) ~ 10.5) simulations and observations are 
nearly consisten t . How ever, one should take into account that 
Kobavash i et al.l d2007l) extract the metallicities in the simulated 
galaxies within a radius of r < 20kpc, i.e. by including metal 
poor external regions, while observations provide metallicities 
within a radius of r < 3kpc. Aperture effects are probably more 
important in large, massive galaxies. If the extraction radius of 



the simulations is matched to the observations then the discrep- 
ancy probably increases strongly also in the high mass region. 

The main problem of the latter simulations seems to be 
that the bulk of the chemical enrichment occurs in small galax- 
ies, yielding a steep metallicity evolution at masses below < 
10 9 M a on the mass-metallicity plane. Then such evolved sub- 
units merge to form massive systems, with little additional star 
formation (close-to-dry merging), implying little additional en- 
richment (i.e. flatter evolution towards high masses). The dis- 
crepancy with the observations seems to imply that the evo- 
lution of galaxies at high redshift occurs through the assem- 
bly of little evolved small galaxies, so that high-z objects with 
log(M*/M Q ) > 9 can exist with relatively low metallicities. In 
these systems most of the star formation and chemical evolution 
occurs once they are already assembled into bigger systems. 

The l atter scenario is well de scribed by the simulations pre- 
sented in lGovernato et al.l (120071) . which model the evolution of 
disks within a hierarchical framework. In these simulations a 
strong feedback, due to SNe and to gas heating by the UV radia- 
tion, prevents small galaxies to evolve significantly before merg- 
ing into a bigger galaxy, while the bulk of the chemical evolution 
and star formation occurs in the gas wh ich has already settled 
into the proto-disk. iBrooks et al.l d2007l) inferred the evolution 
of the mass-metallicity relation from these simulations, whose 
prediction at z=3 is shown in the bottom-left panel of Fig. [TO] 
(Brooks, priv. comm.). The agreement with the observations is 
good, although at high masses the simulations tend to under- 
reproduce t he observed metallic ity. However the metallicities 
provided by Brooks et al] (120071) are obtained without any con- 
straints on the aperture (i.e. all cold gas); if the information from 
the simulated galaxies is extracted within our observational aper- 
ture (r < 3 kpc) then metallicities are expected to increase (espe- 
cially in large massive galaxies) and bet ter reproduce the obser- 
vatio ns. Yet, a po t ential problem of the iGovernato et all ([2007 ) 
and IBrooks et al.l (|2007) models is that their forming disks are 
characterized by modest star formation rates, never exceeding 
~ 20 M yr~\ while LBGs (except for a minority of them) are 
characterized by significantly larger SFR, suggesting that they 
are in the process of rapidly forming spheroids. 

Detailed predic tions on the mas s -meta llicity relation were 
also obtained by Finl ator & Pavel (120081) who used three- 
dimensional hierarchical simulations along with detailed out- 
flows models. They show that the evolution of the mass- 
metallicity relation out to z~2 can be well reproduced if a 
"momentum-driven wind" model is incorporated. The predic- 
tions of their model at z~3 reproduce reasonably well also the 
metallicity in massive galaxies (log(M,/M ) ~ 10.5) observed 
by us. However, their model predicts an up-turn of the mass- 
metallicity slope at z~3 which is not observed by us. Yet, the 
slope of the mass-metallicity relation is still poorly determined 
in our data, due to the shortage of low-mass galaxies; we should 
wait for the completion of the AMAZE program before claiming 
any significant inconsistency with the model in terms of slope of 
the relation. 

In the bottom-right diagram of Fig. [10] we also show the 
mass-metallicity relation expected by the double-infall mod- 
els fo^_fhe_Jorrrmtion_of gala ctic disks and dwarf s presented 
in IChiappini et all d2001l) and ICescutti et all d2007l) . The green 
points show the mass-metallicity relation of spirals according 
to these models, by tracing back their evolution until z=3.5. 
The figure shows a significant discrepancy between the model 
and the observations. Such a discrepancy is not surprising, since 
LBG's at z>2 are probably spheroids in the process of rapidly 
forming stars, and not spirals. More specifically, in the double- 
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Fig. 11. Evolutionary tracks of individual galaxies on the mass- 
metallicity diagram accordin g to some models fo r the forma- 
tion of sphe roids presented in Granato et al] d2004l) (green solid 
lines) and in lPipino et al.l (l2006h (violet dashed lines), compared 
with the mass-metallicity relation observed at z~3.5 (diamonds 
with errorbars). Squares mark the models temporal evolution 
from 0.02 Gyr to 0.3 Gyr, in steps of 0.04 Gyr. Red symbols 
show the masses of the observed objects inferred with the BC03 
templates, while orange symbols are the masses inferred with 
the M05 templates. 



infall model for disks the star formation rate never exceeds a few 
times 1 M yr~ 1 , while the median SFR of LBGs is 90 M Q yr~ 1 . 

Summarizing, there are currently no models or simulations 
that can satisfactory explain the mass-metallicity relation ob- 
served at z~3. The closes t matc h is probably with the simu- 
lations of iGovernato et al] (120071) and iBrooks et al] d2007l) . al- 
though even in these cases there are some discrepancies in terms 
of SFR. The location of LBG's at z~3 on the mass-metallicity 
plane, along with comparison with these models, suggest that 
z~3 galaxies have been assembled through low mass systems 
whose star formation efficiency was suppressed, hence which 
were little evolved. The bulk of the star formation and of the 
chemical enrichment occurred once small galaxies were already 
assembled into bigger systems. In other words, most of the merg- 
ing occurred before most of the star formation. 

We conclude this section by comparing in Fig. [TT] the loca- 
tion of z~3.5 galaxies on the mass-metallicity diagram with the 
evolutionary tracks (as a function of time) expect ed for individ- 
ual sph eroi dal galaxies, accord ing to the models in Grana toet al] 
(2004) and Pip ino et al] d2006l) . These models prescribe a nearly 
monolithic formation of elliptical galaxies, where pristine gas 
collapses from the halo. In these models star formation is grad- 
ually quenched as the galaxy evolves due to the feedback in- 
troduced by star formation and/or AGN activity. The mass- 
metallicity tracks of these models are shown in Fig.QT]for differ- 
ent final st ellar masses: green s olid lines and viole t dashed line 
are for the iGranato et al] (120041) model and for the Pipin o~eTaT] 
(2006) model, respectively. The observational data from the 
Amaze project are shown with diamonds and errorbars. Red 



points are for masses inferred by using the BC03 templates, or- 
ange points are for masses inferred by using the M05 templates. 
The models can easily embrace the observed data points. This 
comparison is not aimed at explaining the mass-metallicity re- 
lation at z~3.5, since models of individual galaxies at different 
masses must be convolved with the evolution of cosmic struc- 
tures to obtain a prediction of the mass-metallicity relation at any 
epoch. Nonetheless, Fig.fTTIshows that the combination of mass 
and metallicities observed in individual star forming galaxies at 
z~3.5 does not lie in a region difficult to populate by models of 
individual galaxies; on the contrary, individual observations can 
be easily explained in terms of rapidly evolving, massive sys- 
tems, through these simple models. However, convolving these 
models with the hierarchical growth of dark matter structures 
is required to verify whether they can really explain the mass- 
metallicity relation at z~3.5. 

8. Summary and conclusions 

We have presented initial results of the AMAZE project, an ESO 
large programme aimed at determining the mass-metallicity re- 
lation of star forming galaxies at z>3. Near-IR spectra are being 
obtained with SINFONI, the VLT near-IR integral field spec- 
trometer, for a sample of Lyman Break Galaxies at 3<z<5. 

Gas metallicities are inferred by using a combination of diag- 
nostics involving nebular lines observable in the H and K bands. 
To have a metallicity scale consistent with the results obtained 
by previous surveys at lower redshifts, we derived new accurate 
calibrations of various strong-line metallicity diagnostics span- 
ning the wide range 7.1 < 12 + log(0/H) < 9.1. 

AGNs (which would affect and make unusable the metallic- 
ity diagnostics) were carefully removed through a multiwave- 
length approach using X-ray, optical and mid-IR data. In partic- 
ular, the use of mid-IR (Spitzer-MIPS) data allows us to discard 
even heavily obscured, Compton thick AGNs. 

Stellar masses are inferred by fitting multi-band photomet- 
ric data with galaxy templates. Within this context crucial is the 
use of Spitzer-IRAC data, which sample the rest-frame near-IR 
stellar light at 3<z<5. 

In this paper we have presented results from an initial sample 
of 9 LBGs at 3.1<z<3.7. Emission lines required to constrain 
the gas metallicity are detected in all sources. From the mass- 
metallicity relation at z~3.5 inferred from this initial sample we 
obtain the following results: 

1 . When compared with the results of previous surveys at lower 
redshift (recalibrated to match our metallicity and mass cal- 
ibrations), we obtain evidence for a clear evolution of the 
mass-metallicity relation as a function of redshift. Since at 
different redshifts we are sampling different populations of 
galaxies, the observed evolution should not be considered the 
evolutionary pattern of individual galaxies, but as the evolu- 
tion of the mass-metallicity relation of the dominant popu- 
lation of star forming galaxies at each epoch. The evolution 
of the mass-metallicity relation is faster at 2.2<z<3.5 than 
at later epochs (z<2.2). This result indicates that z~3 is an 
epoch of major action for the enrichment of galaxies, even in 
high mass systems 

2. There are indications that the metallicity evolution is not 
constant with mass: low mass galaxies evolve more strongly 
than massive systems. If confirmed, this result can be consid- 
ered as the chemical version of the galaxy downsizing, i.e. a 
scenario where high mass galaxies reach high metallicities 
at high redshift, on short timescales, while low mass systems 
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enrich their ISM over a prolonged period of time, extending 
to the current epoch. However, this finding is subject to un- 
certainties due to low statistics at low masses, and must await 
additional data to be confirmed. 
3. The observed mass-metallicity relation at z~3 is difficult to 
reconcile with the simulations by some hierarchical models, 
which predict metallicities higher than observed. If aperture 
effects are taken into account then the discrepancy is proba- 
bly even higher. The main problem of these models seems to 
be that galaxies are assembled once they are already evolved, 
from the stellar and chemical point of view. Our results sug- 
gest that galaxies at z>3 are assembled mostly with relatively 
un-evolved small galaxies, whose star formation efficiency is 
low. Most of the chemical evolution (hence most of the star 
formation) must occurs once small galaxies are already as- 
sembled into bigger ones. This implies that most of the merg- 
ing occurs before most of the star formation. Indeed, models 
and simulations where a strong feedback keeps star forma- 
tion low in the assembling small galaxies (preventing them to 
evolve strongly from the stellar and chemical point of view), 
provide a much better description of the mass-metallicity re- 
lation at z~3.5. 

The finding that galaxies at z>3 are mostly assembled with 
un-evolved sub-units is not necessarily in contrast with mod- 
els of "dry-merging", i.e. models where galaxy assembly oc- 
curs through systems that are significantly evolved and with 
little residual gas. Indeed, dry merging may be the main 
mode of galaxy evolution at lower redshifts (z<3) without 
being in conflict with our findings at z>3. 
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Table 1. Photometric properties and integration times of the AMAZE sub-sample presented here 



Name 


RA(J2000) 


Dec(J2000) 


z 


Rab 


[3.6yum]AB 


[24yum] A B 


Tjnt 


CDFa-C9 


00 


53 13.7 


+ 12 32 11.1 


3.2119 


23.99 


22.66 


>19.3 


3.3 


CDFS-4414 


03 


32 23.2 


-27 51 57.9 


3.4714 


24.18 


22.99 


20.7 


6.7 


CDFS-4417 


03 


32 23.3 


-27 51 56.8 


3.4733 


23.42 


22.23 


>19.9 


6.7 


CDFS-6664 


03 


32 33.3 


-27 50 07.4 


3.7967 


24.80 


25.26 


>19.9 


5.0 


CDFS- 16767 


03 


32 35.9 


-27 41 49.9 


3.6241 


24.13 


24.04 


>19.9 


7.5 


CDFS-11991 


03 


32 42.4 


-27 45 51.6 


3.6110 


24.23 


23.91 


>19.9 


7.5 


CDFS-2528 


03 


32 45.5 


-27 53 33.3 


3.6877 


24.64 


24.17 


>19.9 


4.2 


SSA22a-M38 


22 


17 17.7 


+00 19 00.7 


3.2940 


24.11 


21.75" 


>19.6 


6.7 


SSA22a-augl6M16 


22 


17 30.8 


+00 13 10.7 


3.2920 


23.83 


22.87 


>19.8 


4.2 



The following quantities are reported in each column: col. 1, object name; cols. 2-3, coordinates (J2000); col. 4, redshift spectroscopically deter- 
mined through the [OIII]5007 line in our spectra; col. 5, R-band AB-magnitude; col. 6, Spitzer-IRAC 3.6/jm AB magnitude; col. 7, Spitzer-MIPS 
24yum AB magnitude (lower limits are at 3<x); col. 8, on-source integration time (in unit of hours). 
Note: " This is the AB magnitude at 4.5/mi, since photometry at 3.6/jm is not available. 

Table 2. Physical properties of the sample inferred from their SED 



Name logM, (BC03) logM, (M05) M B U-V SFR E(B-V)„ age 

[Me] |M I [mag] [mag] [M yr'] [mag] [Gyr] 

CDFa-C9 10.18!"-* 10.13™ : 2236 ^024 265 - 25+ oiq °- 06 -nm 

CDFS-4414 10-57!™ 10.48+ggl -22.47 -0.02 113 0.20!™ 0.39!™ 

CDFS-4417 !0-29!° :3 J 10.32!^ -22.79 -0.35 438 - 25 -o?o 0m+ -fn 

CDFS-6664 9.49!^° 9 - 30 -oos " 20 - 93 "°- 70 35 al0 -QQ4 al0 -QQ5 

CDFS-16767 10-05!™ 9 - 90 -o'i3 - 2L9 ° "°- 49 84 al5 -gi al2 S 

CDFS-11991 9 - 69 -(U5 9M qII " 2L64 "°- 46 55 0A0 -lfi 0A0 -m 

CDFS-2528 9 - 76 S 9 - 71 -ooo " 21 - 57 "°- 33 101 0- 20 -oo5 - 06 -ooi 

SSA22a-M38 H.01+™ 10.71+™ -22.99 0.22 115 0.20+™ 0.63+° :95 



0.41 '"'"4« w.^.i/_ 020 

+0.20 in i =+0.11 n An nnA--""i 



SSA22a-augl6M16 10.29!™' jj^j!°^ -22.32 -0.23 42 0.06!^-^ 0.39_ n: , ; 

The following quantities are reported in each column: col. 1, object name; col. 2, stellar mass inferred by using the galaxy templates of BC03; col. 
3, stellar mass inferred by using the galaxy templates of M05; col. 4, absolute B magnitude; col. 5, rest- frame U-V color; co l. 6, star formation 
rate (by using BC03); col. 7, dust reddening affecting the stellar light, by using the attenuation curve o f ICalzetti et all ([2000); col. 8, age of the 
stellar population. 



Table 3. Line fluxes and metallicities inferred from the near-IR spectra 



Name 


F([OIII]5007) 


F(Hy3) 

10-' 7 erf 


F([OII]3727) 

; s -1 cm" 2 ' 


F([NeIII]3870) 


12+log(0/H) 


CDFa-C9 


6.83+0.21 


1.44+0.18 


2.16+0.14 


0.53+0.13 


8.10!™ 


CDFS-4414 


0.95+0.11 


0.40+0.10 


0.87+0.08 


<0.20 


8.54+™ 


CDFS-4417 


2.00+0.11 


0.89+0.09 


1.17+0.11 


0.20+0.09 


o cc+8:oii 

-Q.1Q 
7 711+O.3S 
'• -0.29 


CDFS-6664 


2.71+0.20 


0.50+0.16 


0.40+0.07 




CDFS-16767 


1.90+0.14 


0.50+0.10 


0.50+0.06 




Q Q 1+0.11 

8 -0.42 
8.02!;' 


CDFS-11991 


2.94+0.11 


0.26+0.10 


0.62+0.08 


0.15+0.08 


CDFS-2528 


1.87+0.26 


0.32+0.18 


0.67+0.16 




8.07+1 


SSA22a-M38 


5.41+0.31 


1.48+0.20 


3.51+0.37 


<0.70 


8.34!^ 


SSA22a-augl6M16 


1.65+0.10 


0.30+0.07 


0.45+0.11 


<0.25 


7 (jcj+0.26 


Composite" 


4.58+0.12 


1.00+0.10 


1.50+0.14 


0.25+0.12 


8 24+^ : '^ 



The following quantities are reported in each column: col. 1, object name; col. 2-5, emission line fluxes; col. 6, gas metallicity. 

Notes: " In the case of the composite spectrum fluxes are normalized to the H/3 flux (which is also subject to an error as listed in the corresponding 

column). 



Table 4. Coefficients for different strong-line metallicity diagnostics in Eq.Q] 



Flux ratio (log R) 


t'o 


Cl 


t'2 


C3 


c 4 


logR 23 a 


0.7462 


-0.7149 


-0.9401 


-0.6154 


-0.2524 


log [F([NII]6584)/F(Ha)] 


-0.7732 


1.2357 


-0.2811 


-0.7201 


-0.3330 


log [F([OIIl]5007)/F(H/3)] 


0.1549 


-1.5031 


-0.9790 


-0.0297 




log [F([OIII]5007)/F([OII]3727)] 


-0.2839 


-1.3881 


-0.3172 






log [F([OII]3727)/F(H/J)] 


0.5603 


0.0450 


-1.8017 


-1.8434 


-0.6549 


log [F([OIII]5007)/F([NII]6584)] 


0.4520 


-2.6096 


-0.7170 


0.1347 




log [F( [Nelll] 3870) /F( [Oil] 3727)] 


-1.2608 


-1.0861 


-0.1470 







Notes: " R 23 = (F([OII]3727) + F([OIII]4959) + F([OIII]5007))/F(H/J) 
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Table 5. Best fit parameters for analytical form of the mass-metallicity relation in Eq.|2]at different redshifts 



z 


logMo 


Ko 


0.07 


11.18 


9.04 


0.7 


11.57 


9.04 


2.2 


12.38 


8.99 


3.5" 


12.76 


8.79 


3.5* 


12.87 


8.90 



Notes: " Value obtained by using the masses estimated with the BC03 templates. * Value obtained by using the masses estimated with the M05 
templates. 



